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ABSTRACT 

We investigate the nature of OGLE05155332-6925581, one of the brightest members of 
the enigmatic group of Double Periodic Variables (DPVs) recently found in the Magellanic 
Clouds. The modeling of archival orbital light curves (LCs), along with the analysis of the 
radial velocities suggest that this object is a semi-detached binary with the less massive star 
transferring matter to the more massive and less evolved star, in an Algol-like configuration. 
We find evidence for additional orbital variability and Ha emission, likely caused by an ac- 
cretion disc around the primary star As in the case of (3 Lyr the circumprimary disc seems to 
be more luminous than the primary, but we do not detect orbital period changes. We find that 
the LC follows a loop in the color-magnitude diagram during the long cycle; the system is 
redder when brighter and the rising phase is bluer than during decline. Infrared excess is also 
present. The source of the long-term periodicity is not eclipsed, indicating its circumbinary 
origin. Strong asymmetries, discrete absorption components (DACs) and a 7 shift are new 
and essential observational properties in the infrared HI lines. The DACs strength and RV 
follow a saw-teeth pattern during the orbital cycle. We suggest that the system experiences 
supercycles of mass outflow feeding a circumbinary disc. Mass exchange and mass loss could 
produce comparable but opposite effects in the orbital period on a long time scale, resulting 
in a quasi-constancy of this parameter. 

Key words: binaries: eclipsing, stars: early-type, stars: evolution, stars: mass-loss, stars: 
variables-others 



1 INTRODUCTION: ABOUT DPVS AND 
OGLE05155332-6925581 

The star OGLE05155332-6925581 (MACHO IDs 79.5739.5807 
and 78.5739.78; OGLE LMC-SC8-125836) is a member of the 
group of close binaries named Double Periodic Variables (Men- 
nickent et al. 2003, Mennickent et al. 2005a). These stars show two 
closely related photometric periodicities; the shorter one reflects 
the binary period but the cause for the longer periodicity (or quasi- 
periodicity since it is not exactly cyclic, Mennickent et al. 2005b) 
still is unknown. We have selected OGLE05155332-6925581 for 
a detailed monitoring and study since it is relatively bright, has 
a large amplitude of the long term changes and its eclipsing na- 
ture could help to shed light on the nature of the phenomenon of 



* E-mail: rmennick@astro-udec.cl. Based on observations canied out at 
ESO and CTIO telescopes: ESO proposal 076.D-0126(A) and NOAO pro- 
posal CHILE03B-01 



Double Periodic Variables, especially their evolutionary stage and 
long-term variability. 

OGLE05155332-6925581 shows a long term periodicity of 
188 day and an eclipsing variability with orbital period 7.2843 days 
(Mennickent et al. 2003). The standard multi-color photometry of 
this system is available in several catalogues published recently. 
Zaritsky et al. (2004) gives simultaneous U = 15.040 ± 0.032, B 
= 15.575 ± 0.031 and V = 15.587 ± 0.057 mag. Wyrzykowski 
et al. (2003) gives V = 15.38, B = 15.32 and I = 15.27 mag 
at maximum. Ita et al (2004) provides single-epoch simultane- 
ous three color infrared imaging with J = 15.586, H = 15.566 
and Ks = 15.490 mag. No epoch is available for these measure- 
ments. OGLE05 155332-6925581 is one of the brightest DPVs in 
the LAIC, in consequence the quality of the available photometry 
is good enough for concluding about physical properties of this ob- 
ject. On the other hand, their photometric credentials are typical for 
DPVs so it can be considered as a representative case. 

OGLE05 155332-6925581 appears in the field as a visual 
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Figure 1. The finding chart of OGLE05155332-6925581 based on an I- 
band OGLE image. North is up and east right. The field of view is 1 f 7 X 1 f 7 
and centered at RA(2000) = 05:15:53.32 and DEC(2000) = -69:25:58.1. 
The circle indicates the position of OGLE05155332-6925581. 



double, with a nearby companion 3.8"to the East. This neighboring 
star was also observed photometrically and spectroscopically. The 
OGLE position for the companion is RA = 05:15:53.88 and DEC 
= -69:25:58.1 and the photometry gives V = 16.947, (B - V) 
= 1.812 and {V - I) = 1.821 mag. The 7-band light curve (LC) 
consisting of 395 datapoints gives a mean of / = 15.072 mag with 
RMS 0.006 mag. 2MASS data for this star give J = 13.928 ± 
0.048, H = 13.046 ± 0.061 and K = 12.817 ± 0.055 mag. A find- 
ing chart based on an OGLE image taken at the /-band is given in 
Fig.l. 

We provide a description of our observations and archive pho- 
tometric data in Section 2 along with our methods of data reduction 
and analysis. The results of our LC modeling and spectroscopic 
study are given in Section 3 and a discussion is presented in Sec- 
tion 4. Our conclusions are given in Section 5. 

In this paper we will label with subindex 1 those parameters of 
the more massive star (the primary) and with 2 those of less massive 
star of the binary pair. 



2 SPECTROSCOPIC OBSERVATIONS AND SOURCES 
OF PHOTOMETRY 

We have monitored OGLE05155332-6925581 spectroscopically 
in the optical wavelength range during August 2003 and February 
2004, and in the infrared region during December 2005 and De- 
cember 2006 and also obtained an optical high resolution spectrum 
at a single epoch in 2007. In this section we give details of these 
observations. In addition, we have used public domain photomet- 
ric data available from the MACHO (Allsman & Axelrod 2001) 
and OGLE II databases (Udalski et al. 1997, Szymanski 2005) 
along with OGLE III data. The MACHO light curves analyzed in 
this paper consist of 1469 photometric measurements obtained at 
the "red" filter (central wavelength 7000 A) and 1555 measure- 
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Figure 2. The average infrared spectra for the nearby companion observed 
in the same slit. Equivalent total exposure time is 3.3 hours. The spectra 
have been coirected by the stellar radial velocity. Line identification is based 
on the work by Wallace et al. (2000). 



ments obtained at the "blue" filter (central wavelength 5200 A) dur- 
ing HJD 2 448 826. 136 and 2451 536.193 (red) and 2451 541.999 
(blue). We also analyzed 896 OGLE II and III magnitudes obtained 
at the I-band (central wavelength 8000 A) spanning a time interval 
between HJD 2450457.649 and 2453 823.498. The total time cov- 
ered by the photometric data analyzed in this paper is 13.7 years. 

2.1 Time resolved optical spectroscopy 

26 spectra of OGLE05 155332-6925581 in the wavelength re- 
gion of 3900—5300 A were obtained with the 1.5m telescope 
of the Cerro Tololo Interamerican Observatory (CTIO), Chile, in 
service mode using the Cassegrain Spectrograph with grating 26 
tilted 15?94. The slit width of l"yielded a nominal resolution 
of 2 A. He-Ar comparison spectra obtained along every science 
exposure provided wavelength calibration functions with typical 
RMS of 0.1 A (6 km s~^ at H/3). Standard techniques of spec- 
trum extraction, calibration and analysis were applied using IRARj 
tasks. We measured the radial velocity (RV) of the Balmer lines 
with simple gaussian fits. The average RV of the standard star 
LTT2415, measured every observing night, resulted to be 255 ± 
36 km s~^, which compares well with the value of 253.4 km s~^ 
reported in the SIMBAD database. Our measured velocities for 
OGLE05 155332-6925581 are given in Table 1. 



2.2 The MIKE high resolution spectrum 

We obtained one spectrum of OGLE05155332-6925581 with the 
Magellan Inamori Kyocera Echelle (MIKE) in the Baade Magel- 
lan telescope of Las Campanas Observatory, Chile, the night of 
February 23 (UT), 2003 (UT-start 00:58:39, HJD at mid-exposure 
2452 693.55208). We used a 1800 sec exposure and a slit width 



^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 3. Results of the Fourier light curve decomposition showing the orbital (upper graphs) and long-term (lower graphs) variability separately. 



Table 1. The radial velocities for OGLE05155332-6925581 in km s^^ 
referred to the Local Standard of Rest. HJD' means HJD-2 450000. The 
orbital and superphases are given along with detection d of He 1 4471 (p: 
positive, w: weak, w: weak doubtful and n: negative). 



UT-date 


HJD' 


*o 


$„ 


H(5 


H7 


H/3 


d 


Aug08 


2859.8912 


0.582 


0.074 


226 


154 


208 


n 


Aug20 


2871.8437 


0.222 


0.144 


- 


- 


- 


n 


Sep03 


2885.7998 


0.138 


0.225 


293 


308 


328 


n 


Sep 11 


2893.7845 


0.235 


0.271 


366 


410 


357 


n 


Sep24 


2906.8188 


0.024 


0.346 


359 


366 


312 


w 


Sep26 


2908.8315 


0.300 


0.358 


508 


295 


459 


P 


OctlO 


2922.8309 


0.222 


0.439 


408 


430 


316 


w 


Oct22 


2934.7777 


0.862 


0.508 


207 


212 


185 


p 


Oct24 


2936.7097 


0.127 


0.520 


142: 


320 


240 


w 


Nov27 


2970.8199 


0.810 


0.717 


132 


229 


215 


w: 


Nov28 


2971.7277 


0.935 


0.723 


225 


191 


231 


P 


Dec 10 


2983.7592 


0.586 


0.792 


101 


136 


244 


P 


Dec 11 


2984.7069 


0.716 


0.798 


205 


103 


241 


n 


Dec 12 


2985.6706 


0.849 


0.803 


219 


164 


165 


n 


Dec 17 


2990.7294 


0.543 


0.833 


217 


306 


274 


n 


Dec 18 


2991.7520 


0.684 


0.839 


235 


269 


148 


w: 


Decl9 


2992.7313 


0.818 


0.844 


109 


166 


155 


w: 


Dec30 


3003.6639 


0.319 


0.908 


346 


378 


397 


n 


Dec31 


3004.6364 


0.452 


0.913 


290 


274 


248 


w: 


Janll 


3015.6778 


0.968 


0.977 


267 


268 


241 


n 


Janl2 


3016.6215 


0.098 


0.983 


- 


467 


448 


n 


Jan20 


3024.7157 


0.209 


0.030 


391 


451 


471 


n 


Jan21 


3025.6757 


0.341 


0.035 


391 


267 


405 


P 


Jan22 


3026.6209 


0.471 


0.041 


256 


253 


347 


P 


Jan31 


3035.6290 


0.707 


0.093 


116 


182 


95 


P 


FebOl 


3036.6371 


0.846 


0.099 


164 


67 


131 


P 



of 0.35", resulting in a resolving power ~ 70000. Our spectrum 
was taken at orbital phase 0.756 and superphase 0.110 (see below). 
The calibration and extraction of the spectrum were performed with 
standard IRAF tasks. The spectrum was wavelength calibrated with 
several hundreds comparison lines of ThAr spectra taken before 



and after the exposure. The RMS of the calibration function was 
about 10"^ A. 



2.3 J-band infrared spectroscopy with the ESO ISAAC 

Our observations were conducted at the European Southern Obser- 
vatory (ESO) Paranal observatory, Chile, with the Infrared Spec- 
trometer and Array Camera (ISAAC) in the SW-MR mode in 
service mode (http://www.eso.0rg/instruments/isaac/l. The detec- 
tor was the Hawaii Rockwell array with nominal gain 4.5 e/ADU 
and readout noise II electrons. This spectroscopic mode provided 
a plate scale of 0. I47"per pixel. We used a slit width of 1 "yielding 
a spectral resolution of 2.4 A and resolving power R ~ 4550. We 
observed OGLE05 155332-6925581 at three different spectral re- 
gions: thej region between 1.061 — 1.108 /xm, the sz region between 
1.154—1.215 /im and the jplus region between 1.265—1.324 /im. 
The observations were done as usual in the infrared, with the tech- 
nique of dithering. This method consists in observing the object 
at two or more positions along the slit. The sky is effectively re- 
moved by subtracting one frame from the other, registering the two 
beams and then subtracting again. The first steps of image reduc- 
tions like geometrical distortion correction, flat fielding and bias 
removal were done by the ISAAC pipeline software at the obser- 
vatory. The wavelength calibration was done identifying night sky 
lines in the image and building a wavelength calibration function 
for every spectrum. Extraction was done with the usual IRAF tasks. 
We used a generic template of infrared atmospheric lines to re- 
move telluric lines. This template was rebinned and resolution de- 
graded to match the characteristics of our spectra. Then we divided 
the science spectra by the telluric spectrum after finding the best 
wavelength shift and vertical scaling providing the best telluric fea- 
ture removal. Equivalent widths of the absorption lines were mea- 
sured by line fitting with Lorentz functions. For those lines showing 
discrete absorptions in their wings, we deblended the profile with 
two Lorentz functions of variable depth and width to take into ac- 
count the blending. For comparison, we also measured the equiva- 
lent width with the standard formula integrating the line flux under 
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Table 2. Observing log for the infrared spectroscopy. We give the exposure time 
in seconds At, the observing block name, the seeing value as measured in the 
acquisition image (si) and read it from the image header (s2), in arcseconds, and 
the spectral band. NA means not available. 



Table 5. Equivalent widths and radial velocities in the infrared spectra. $<, 
is the orbital phase. 



OB 



<I>o EW(A) RV(kms-i) 



line 



date(UT) 


At/4 


MJD 


OB 


si- 


-s2 


band 


21-22/12/05 


600 


3726.2744-3726.2958 


051221 


0.64- 


-0.93 


J 


10-11/01/06 


600 


3746.1699-3746.1913 


205091 


0.54- 


-0.68 


J 


10-11/01/06 


600 


3746.1992-3746.2205 


205094 


NA- 


-0.69 


sz 


11-12/01/06 


500 


3747.1895-3747.2074 


205088 


0.58- 


-0.79 


jplus 


11-12/01/06 


600 


3747.2141-3747.2355 


205092 


NA- 


-0.56 


j 


15-16/02/06 


600 


3782.1062-3782.1276 


205095 


0.60- 


-0.69 


sz 


16-17/02/06 


600 


3783.0996-3783.1209 


205097 


0.57- 


-0.92 


jplus 


14-15/04/06 


600 


3839.9825-3840.0038 


205093 


0.69- 


-0.78 


j 


17-18/03/06 


500 


3811.9970-3812.0150 


205090 


0.60- 


-0.83 


j 


12-13/10/06 


600 


4021.2740-4021.2954 


205098 


0.75- 


-1.32 


jplus 


09-10/11/06 


600 


4049.3222-4049.3435 


205099 


0.56- 


-0.45 


jplus 


06-07/12/06 


600 


4076.2686-4076.2899 


255548 


0.66- 


-1.36 


jplus 



051221 


0.599 


3.5 


205091 


0.331 


5.1 


205094 


0.336 


- 


205088 


- 


5.5 


205092 


0.475 


5.7 


205095 


0.264 


- 


205097 


0.401 


7.1 


205093 


0.210 


4.2 


205090 


0.368 


4.9 


205098 


0.098 


5.8 


205099 


0.948 


5.9 


255548 


0.648 


4.6 



330 

535 



334 



Pa7 
Pa7 



Pa7 



447 


Pa/3 


559 


Pa7 


467 


Pa7 


461 


Pa/3 


279 


Pa/3 


229 


Pa/3 



Table 3. Periods and amplitudes found in the light curves. 



period (d) range (mag) MACHO period (d) range (mag) OGLE 



7.28428 


0.773 


blue 


7.28425 


0.660 


11 


188.4 


0.166 


blue 


185.6 


0.240 


II 


7.28431 


0.713 


red 


7.28430 


0.651 


111 


188.3 


0.206 


red 


172.5 


0.220 


111 



the interpolated continuum, yielding results about 15% lower than 
those obtained with the first method. For the OB205093 spectrum, 
only this method was possible due to the large line asymmetry. We 
measured radial velocities in the infrared spectra determining the 
wavelength at minimum of the flux. We determined an internal er- 
ror of 10 km s^^ from radial velocities of the comparison star, red 
companion described in Section 1, placed in the same slit. Details 
of the observations are given in Table 2. We found that the com- 
parison star is a probable member of the LMC, its RV referred to 
the LSR is -1-246 ±10 km s~^ . The spectra for this star are plotted 
in Fig. 2, along with line identifications; the star could be consid- 
ered as a flux and radial velocity standard in the LMC. Using the 
ratio between the Mn I feature strength to that of Mg I in the for- 
mula provided by Wallace (2000) for Galactic giant stars we obtain 
a temperature of 3953 K. We also note that the Hipparcos ISAAC 
telluric standards Hip 026368 and Hip 029635 are Be stars; they 
show double peak Pa/3 emission with peak separation of 205 and 
235 km s^^, respectively. The first star is the well known Be star 
AZDor (B9.5 Ve) and the second one (HD 44533, B8 V) is a new 
Be star. 



3 RESULTS 

3.1 Photometry from OGLE and MACHO 

The fight curves were analyzed with a Fourier disentangling 
method. This allowed us the separate the Fourier amplitudes 
and frequencies of the two different periodic signals detected in 
the time series photometry. A summary of these parameters for 
OGLE05 155332-6925581 is given in Table 3. The results of our 
analysis of the OGLE and MACHO photometry are displayed 
in Fig. 3. The light curves are folded with the orbital and the 
long-term period. The double eclipsing nature of the variable is 



evident as well as the hump-like character of the long-term light 
curve, whose variability is limited to half a cycle only, and shows 
larger amplitude at the red band. We observe the eclipses slightly 
deeper in the blue band. The improved ephemeris for the main 
minimum is (the numbers in parentheses are the errors with the 
zeros omitted): 

HJDo = 2 450 000.1392(21) + 7.284297(10) x E. (1) 

We note that the long-term variability source is not eclipsed. 
This result contradicts the finding by Mennickent et al. (2005a, 
hereafter M05) based on the visual inspection of the light curves 
phased with the orbital period for two DPVs (their Fig. 12). The 
apparent absence of long-term variability during the main eclipse 
turned to be only a visual effect produced by the large magnitude 
gradients found around minimum. The better treatment of the data 
by the Fourier decomposition technique completely separates the 
short and long-term variability and reveals that the visual impres- 
sion is not real. Our finding suggests that the source of the long- 
term variability is placed around the binary system, and this is the 
rule for all eclipsing DPVs (Kolaczkowski &. Mennickent in prepa- 
ration). 

In Fig. 4 we show the O-C diagram for times of maximum 
of the long-term variability. The maximum moments and their 
errors were measured at every cycle by fitting the LC with the 
average light-curve defined for a given band as a spline function. 
There is a clear change in behavior after cycle number six; the long 
variability becomes faster, and the O-C diagram is consistent with 
a period change from 188.2 d — > 172.5 d. Our data does not reveal 
if this transition is gradual or consist of discrete period jumps. 
The shorter period is observed during the OGLE III epoch, and we 
note in Fig. 3 that the shape and amplitude of the long-term light 
curve remains unchanged, besides the notable period change. The 
ephemeris for the maximum of the long-term variability during 
OGLE III epochs is HJDmaa: = 2 452 674.6 + 172.5 x E. 

During minimum and maximum of the long cycle light the 
system has the most blue and red color, respectively. Fig. 5 shows 
the color-magnitude diagram for the long-term variability of the 
MACHO red light curve, indicating magnitudes taken during ris- 
ing or decline from maximum. We observe a clear pattern in the 
diagram: the loop characterized by bluer rising than decline phase. 
This loop is reminiscent of that observed by de Wit et al. (2006) 
in light curves of outbursting Be stars in the Magellanic Clouds, 
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Table 4. The best system parameters (p), derived assuming a stellar origin for the deeper eclipse. At every band we give the flux fraction l\/{l\ + I2) at 
maximum, g is the surface gravity incms"^. We also give the parameters for the best fit RV = '^ + K2Sin{2-K^o)- f is the mass function. For our calculations 
we have used published values of ^ = 18.39 and E{B — V) = 0.17. 



p 


MACHO blue 


OGLE II 


P 


value 


Line 


7(kms-i) 


i^2(kms-i) 


KMq) 


9 


0.287 


0.393 


Spl 


7 


Pa/3 


347(5) 


167(10) 


3.52(63) 


log 31 


3.96 


3.79 


Sp2 


B6III 


Pa7 


380(18) 


152(25) 


2.65(1.30) 


log 92 


2.86 


2.88 


Mv,2 


-2.83 


H/3 


280(12) 


122(16) 


1.37(54) 


Ti{K) 


25 000 


25 000 


^bol,2 


-4.03 


H7 


271(13) 


123(18) 


1.40(61) 


T2{K) 


14090 


14 330 


R2/a 


0.29 


H(5 


277(11) 


134(15) 


1.82(61) 


i(deg) 


80.4 


80.4 


Ri/a 


0.16 










Ib 


0.419 


0.419 


Po{d) 


7.284297(10) 










fv 


0.392 


0.418 


Ps(d) 


188.2-172.5 










fi 


0.360 


0.387 














fj 


0.324 


0.351 














fH 


0.315 


0.342 
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0.335 
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Figure 4. The O-C diagram for the maxima of the long-term variability. 
Filled squares correspond to MACHO-red data, asterisks to MACHO-blue 
data, open squares to OGLE II data and open circles to OGLE III data. 



Figure 5. The color-magnitude diagram for the long-term variability in the 
MACHO red LC. Solid circles are magnitudes taken during the rising and 
open circles are those obtained during decline from maximum. Magnitudes 
taken at the maximum and broad minimum are shown by pluses. The loop- 
like solid line follows the average temporal trend. 



and is also a strong evidence for a circumbinary source for the long 
period variability. 



3.2 Orbital light curve modeling 

We used the Wilson-Devinney (Wilson & Devinney 
1971) code for modeling the eclipsing variability of 
OGLE05 155332-6925581 but we encountered problems 
with getting convergence of the solution for all available LCs. We 
found also non-physical solutions revealed as a negative basal 
constant light. Therefore we decided to use a light curve modeling 
program based on WD code that find the best solution mapping a 
global minimum of summed residuals in a multi parametric space 
using Monte-Carlo techniques (Michalska & Pigulski 2005). 
We conducted the search for an optimal solution for every band 
separately applying also reliable simplifications. Circular orbits, 
synchronous rotation and absence of third light are assumed by this 
method. A fixed temperature for the primary component is roughly 



estimated from the brightness and the distance of the system 
using the method described in Michalska & Pigulski (2005). This 
iterative method picks up the visual magnitude at maximum and 
starts assuming equal luminosities for both components and the 
Ti versus Mtoi relation for Galactic detached main-sequence 
eclipsing binaries published by Harmanec (1988). The difference 
between the Galaxy and LMC metallicities for the parameters BC 
and Mv was calculated from the models of Bertelli et al (1994). 
The Ti value hence estimated converges to 25 000 K. The best 
fit parameters are given in Table 4 and the model light curves are 
shown in Fig. 6. Estimates of the fit quality are given in Figs. 7 and 
8. We observe systematic differences between the best solutions at 
different wavelength bands (Figs. 7 and 8). We argue that these 
differences could be due to the absence of including a constant 
basal light (third light, probably residual light from the source 
of the long-term variability) in our models. In consequence we 
cannot adopt these models as the conclusive physical solution. 
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Figure 6. A detailed comparison of the model and the observations. Red open circles are observations (MACHO blue; long term variability is removed), the 
dark blue line is the fit of spline functions to the observations, the green line is the model, based on the parameters of the best solution for the MACHO blue 
data. Magenta crosses are the observations minus model, the light blue line is the spline fit minus model. Zoomed views of eclipses are given in the right side 
panels. 
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Figure 7. Results of our Monte Carlo simulations. Every point represents a model, characterized by a set of input parameters for the WD code and its sum of 
squares of residuals (ssr). The models in the figures are those minimizing ssr at the MACHO blue LC (upper panels). These models generate LCs at bandpasses 
corresponding to MACHO red and OGLE LCs. Comparison of these synthetic LCs with the observational data yield to the diagrams of the middle and bottom 
panels, respectively. Black dots are a 10% selection of the best models for the MACHO blue LC. 



but only for a discussion concerning the general properties of the 
system. 

The potentials of the secondary star in the whole possible 
range of the mass ratios clearly indicate a semi-detached config- 
uration (Fig. 9). We find Ti = 25 OOOiC, T2 = 14 090i^, log gi = 
3.96 {g is the surface gravity in cm s~^), \0gg2 ~ 2.86, i — 



80? 4 and q = 0.29 from the best MACHO blue light curve fitting 
model. These parameters are well confirmed by the analysis of the 
OGLE II light curve except the mass ratio, for which we find q = 
0.39. The ratio of the stellar radii R2/R1 is in the range 1.79 (for 
OGLE Il-band data) to 1.91 (for MACHO blue). 

The temperature and surface gravity of the secondary star 
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Figure 8. Same as Fig. 7 but considering models optimized for tlie OGLE II liglit curve. 



agree with the improved values obtained from the spectra reported 
by M05, viz. T2 = 13 500 ± 1000 K, log g2 = 3.2 ± 0.5 and M„ 
= -2.0 ± 0.3, Mtoi = -2.9 ± 0.3 and spectral type B5 IV. 

Further detailed comparison between synthetic and observed 
LCs reveals additional components (fourth light) modulated with 
the orbital cycle of the system. As we can see in Fig. 6 the overall 
model LC is quite good, but some significant deviations are present 
at phases of both eclipses. The system is fainter than the model 
at the main minimum and also during ingress and egress of main 
eclipse. In addition, the observed secondary eclipse is wider than 
the model. These features on the observed LC could be explained 
by the presence of an optically thick disc of gas around the primary 
star. This kind of accretion structure could be present in the sys- 
tem due to mass transfer from the secondary to the primary star. 
The deeper minimum is the opposite effect that expected in the 
presence of the third light. In the observed eclipses both additional 
components, third and fourth light, should have a contribution in 
opposite directions, but in our simple modeling is impossible to es- 
timate their total influence. We will return to the LC interpretation 
in the discussion section. 



3.3 The infrared spectra: DACs and 7 shifts 

We show in Fig. 10 and 11 the smoothed (7-point average) Pa/3 
and Pa7 hydrogen lines for OGLE05155332-6925581. Equivalent 
widths and radial velocities are given in Table 5. The HI lines ap- 
pear in absorption with variable asymmetries and depths; the depth 
is larger at phases 0.37—0.40 and smaller at phases 0.60—0.65. The 
H I profiles looks quite complex in structure. We identify the main 
absorption with the secondary star (see Section 3.4) but there is a 
weaker and wider absorption component in the H I lines appearing 
in some phases which we attribute to the primary star. In Fig. 12 we 
show the non smoothed Pa/3 spectrum at phase 0.47, along with a 
synthetic binary absorption profile built considering the stellar pa- 
rameters found in Section 3.2. In order to reproduce the observed 



Pa/3 profile, we consider the contribution of a primary star with 
Te// = 25 000 K and log g = 4 and a secondary star with Te// = 
14 000 K and log g = 3. For each object we assumed that the photo- 
sphere is in radiative and hydrostatic equilibrium and is represented 
by the corresponding density and temperature distributions given 
by Kurucz's (1979) model atmospheres. The line profiles are ob- 
tained by solving the transfer equation and the equations of statisti- 
cal equilibrium for multi-level atoms in the commoving frame with 
the ETLA code (as described by Catala & Kunasz, 1987). We as- 
sumed a spherically symmetric slowly expanding atmosphere with 
a mass loss rate 10"^ MQ/yr. Once the line source function was 
known, the line profiles for the primary and secondary stars were 
computed in the observer's reference frame and convolved with a 
grid of rotational profiles. Finally, we summed up both contribu- 
tions in order to fit the line profile of Pa/3 of the binary system. It is 
notable the high velocity required to fit the primary star spectrum, 
viz. 400 km s^^. This fact shows that the primary rotates very fast, 
non-synchronously and probably at critical speed. We do not ob- 
serve emission lines in the studied infrared spectral regions. 

We find transient absorption features around A 1.084, 1.093 
and 1.280 fim. We have identified the feature around A 1.084 fim 
with He I 1.083 /im. During orbital cycle this line is quite vari- 
able in strength and shape; it is stronger and wider at phases 0.33 
and 0.37 and disappears in the spectrum of phase 0.21. Their RV 
is confusing, it follows H7 except at phases 0.33 and 0.37 when 
is bluer and stronger. The other two absorption features are much 
narrower and placed at the blue side of the H I lines, with velocities 
between 550 and 850 km s~^ from the line center. The good re- 
moval of telluric lines in the region 1.086—1.092 /im illustrated in 
Fig. 10 suggests that the transient features are real and not artifacts 
of the process of telluric correction. While we are confident about 
the well defined blue-shifted narrow discrete absorption compo- 
nents (DACs), the quality of our data is not enough to be conclu- 
sive about the reality of some narrow emissions features observed 
in some spectra. 



8 Mennickent et al. 



We find DACs stronger on superphases 0.10, 0.76, 0.81 and 
0.97, i.e. around long-cycle maxima, and weaker or absent on su- 
perphases 0.40—0.65, during long-cycle minimum (Fig. 13). This 
finding could indicate that DACs are associated to the long-term 
periodicity. However, the correlation with their RVs is not so sig- 
nificant as with the orbital period (Fig. 13). We observe DACs with 
higher velocities after phases and 0.5, followed by a monotonic 
decrease of their RVs; they trace a saw-teeth curve during every 
orbital cycle. Based on this behavior we interpret DACs as diag- 
nostics of mass outflow from the binary. The pattern observed in 
Fig. 13 can easily be interpreted as DACs arising from gas streams 
emanating from the external Lagrangian points L2 and I/3. If gas 
accelerates through these streams following large open arcs around 
the binary in opposite direction to the binary motion, the observed 
pattern can be reproduced. Just before phase zero, we observe the 
weak (low density) high velocity component of the L-j, stream. Just 
after phase zero we observe the strong (high density) low velocity 
component of the L2 stream. 

We conclude that the absorption lines observed in the infrared 
spectra of OGLE05 155332-6925581 show a complex and vari- 
able structure, completely different to those observed in normal B- 
A type stellar atmospheres (Groh et al. 2007), and that our obser- 
vations can be interpreted as evidence of mass loss in the binary 
system through the L2 and Z/3 points. 

3.4 The optical spectra, radial velocities and stellar masses 

The average SMARTS optical spectrum shows H I and He I absorp- 
tion lines (Fig. 14). We also detect Si 111 4553, Oil 4070 and Oil 
4416. A list of features is given in Table 6. The TiO bands usually 
found in KO— M2 giant/supergiant stars probably arise from the not 
resolved nearby companion included in the slit. The presence of 
O II 4070 lines indicates a forming region of high temperature and 
the line strength ratio He 1 4387/0 II 4416 ~ 1 is indicative of high 
luminosity. This finding suggests the existence of a high tempera- 
ture interacting region as seen in some Algol-type systems. From 
the inspection of the individual spectra we realized that they are 
highly variable: the helium and oxygen lines sometimes are visible 
and sometimes disappear. Although helium lines and O II 4070 ap- 
pear at all superphases, they are almost always detected during su- 
perphases 0.3—0.5 whereas at other superphases they may or may 
not appear. This suggests a dependence of line strength with su- 
percycle for these lines. The equivalent widths of the Balmer lines 
change notoriously; we find that spectra taken around phase 0.5 
show weaker Balmer lines that those observed around phase 0.0, 
consistent with a secondary star responsible for most of the H I ab- 
sorption. No emission lines were detected in the optical region in 
these medium resolution spectra. 

The high resolution MIKE spectrum shows asymmetrical H I 
absorption lines and Ha appears weakened and with a narrow emis- 
sion component of FWHM = 1 .03 A at A = 6568.24 A (Fig. 15). 
It is not clear if the H I lines asymmetry is due to blending of H I 
lines of the primary and secondary star or it is due to the presence 
of a central P-Cygni type absorption-emission minicomponent as 
suggested by the Ha line. We failed to separate components with 
multi-gaussian deblending. No helium line (neither other type apart 
hydrogen) was observed in the MIKE spectrum. 

Radial velocities of the hydrogen lines based on SMARTS 
spectra are displayed in Fig. 16 along with the OGLE photometry. 
To our surprise, we observe that the velocities trace the motion of 
the secondary; the cooler component. This fact is consistent with 
our finding in Section 3.2 of a more luminous secondary star in the 




Figure 9. Mass ratio vs. potential of the secondary component. Only the 
best models are shown (the black points of previous figures). Parameters 
for OGLE fits are shown as solid circles, squares indicates MACHO blue 
models. The black line is a critical potential for object filling Roche lobe. 
Our results of Monte-Carlo simulations indicate a semi-detached configu- 
ration for OGLE05155332-6925581. 



optical. The parameters of the best fitting sinus functions are given 
in Table 4. The 7 velocity for the Balmer lines is ~ 275 km s^^ 
consistent with objects in the LMC. The radial velocities of the 
infrared lines, measured at minimum of absorption profile follow 
the same behavior that the optical lines, but with larger amplitude 
and also with surprisingly large 7. The primary contributes to the 
HI line motion in opposite direction to the secondary star, then K2 
should decrease, and this effect would be larger in low resolution 
spectra, as observed. The fact that the optical and infrared lines dif- 
fers in 7 is puzzling, and cannot be due to instrumental effects, as 
demonstrated by the unshifted position of telluric lines in the in- 
frared raw spectra and the consistent optical RV with two different 
instrumental setups. We note that the 7 shift is quite similar for 
both Paschen lines, cannot be explained by line blending and is sta- 
ble through several months for every line set. We cannot explain 
this 7 shift at present. 

The mass function for a binary in a circular orbit can be 
expressed as: 



K2 



f = nTyl = 1-0361 X lQ-'{-r^^) 

■' q(l + q) ^ kms ^ ' 






(2) 



where K2 is the half-amplitude of the RV of the secondary star and 
Po the orbital period. The / values derived from Eq. (2) and from 
our radial velocity study are listed in Table 4. Considering the range 
of / and q, the highest possible value for M2 is 2.80 ± 0.5 Mq, 
whereas the minimum is 0.7 ± 0.3 Mq. 



3.5 Differential-slit and wide-band photometry 

We performed differential slit photometry aimed to detect the 
amplitude of the long-term variability in the infrared. The ob- 
serving conditions were optimal to compare the spectral flux of 
OGLE05 155332-692558 1 with their nearby companion placed si- 
multaneously in the slit. Our observations were conducted under 
conditions of photometric sky, we used a wide slit and the seeing 
was always below 0.8". We obtained instrumental magnitudes in 
the different spectral regions defined as —2.5 log / where / is the 
sum of spectral counts (ADUs) in the regions j, jplus and sz. The 



A clue for Double Periodic Variables 9 



central wavelengths were chosen at 1.085, 1.295 and 1.184 fim and 
the band widths 0.020, 0.025 and 0.025 ^m. After checking for 
possible slit misalignments and detector noise characteristics we 
realized that the error of these magnitudes is mainly due to photon 
noise, and neglected in all cases. When analyzing the magnitudes, 
we assumed that the short-term variability in the infrared had the 
same amplitude and behavior that in the /-band and removed this 
component obtaining only the long-term infrared magnitudes. Our 
magnitudes are plotted along with the I magnitudes in Fig. 17. A 
vertical shift was applied to minimize the scatter with the / magni- 
tudes. Our data suggest that there is not significant amplitude dif- 
ference between the long variability in the J-band and the /-band. 

It is difficult to make an analysis of the broad-band photom- 
etry due to the composite nature of the spectral energy distribution 
(SED) and the intrinsic variability of OGLE05155332-6925581. 
To have a rough idea about the SED properties we used the UBV 
photometry by Zaritsky et al. (2004) to calculate an effective red- 
dening free parameter Q = —0.526 which gives [B~-V)o = —0.186 
(Johnson & Morgan 1953) and therefore E{B - V) = 0.174. The 
unreddened effective color corresponds to a B4 III star (Fitzgerald 
1970). The color excess is comparable with the expected in the di- 
rection of the LMC, namely E{B - V) =0.17 (Zaritsky 2004). 
On the other hand the infrared colors derived from Ita et al. (2004), 
viz. (J-K) = 0.096 and (H-K) = 0.076, are too red for AB-type 
stars (Koornneef 1983). This result is also revealed in the spectral 
energy distribution through multiwavelength broad-band photom- 
etry (Fig. 18). This figure shows that the optical spectral distribu- 
tion at maximum light is best represented by the 20 000 K model, 
roughly consistent with comparable contributions of both stellar 
components at this phase. The infrared data, although not necessar- 
ily obtained at maximum as the optical one, clearly indicates that 
there is a significant infrared excess. This color excess seems to in- 
crease beyond 2.2 fim, obviously carmot be reproduced by any of 
the binary stellar components and probably comes from the cool 
circumbinary disc responsible of the long-term variability. 

Our models indicate that, at phase zero, the secondary star 
contributes about seven times more light that the primary in the 
MACHO-blue band. If we use the fact that the system has V = 
16.12 at minimum and assume a distance modulus for the LMC 
of 18.39 ± 0.05 (van Leeuwen et al. 2007) and a color excess of 
E{B ~ V) = 0.17, we get My = -2.83 for the secondary star, 
close to the spectroscopic estimate. Among semidetached Algols, 
only V 356 Sgr (B3V-^A1-2III; Budding et al. 2004, Ibanoglu et al. 
2006) has a more luminous secondary. The mass-luminosity rela- 
tionship for main sequence massive stars L ~ M^''^^ (Vitrichenko, 
Nadyozhin and Razinkova 2007) implies in our case a mass of 12.9 
Mq. Since the dynamical mass is much lower, this result suggests 
that the secondary star could be much overluminous and evolved. 

For a star filling their Roche lobe the mean density is con- 
strained by the orbital period for q ^ 0.8 (e.g. Frank, King and 
Raine 2002): 



^nt: 



^2 gem ^ 



(3) 



In the case of OGLE05 155332-6925581 we find p = 3.6 x 10"'^ 
g cm~^ or log p/p^ = —2.59. Approximating R2 with the Roche 
lobe radius of the secondary star (Eggleton 1983): 



0.6g2/3+;„(l+ql/3) 



(4) 



Table 6. Lines detected in the average optical spectrum of 
OGLE05155332-6925581. The symbol f indicates probable lines of 
the nearby companion. 



A (A) 


line 


A (A) 


line 


3908.7 


SiI3906t 


4474.4 


He 1 4471 


3935.2 


Call 3934 


4717.1 


HeI4713 


3973.2 


He 


4537.4 


TiI+TiII4533t 


4105.4 


H5 


4556.0 


Si III 455 3 


4014.5 


TiI+FeI4010t 


4865.5 


H/3 


4028.7 


He 1 4026 


4924.9 


He 1 4920 


4074.6 


on 4070 


4969.1 


MnI4966t 


4146.9 


HeI4143 


5003 


TiO band 5003t 


4344.4 


H7 


5170 


TiOband5167t 


4390.4 


He 1 4388 


5274.8 


CaI+FeI5270t 


4420.3 


OII4416 








we found R2 = 0.28a or 0.30a, for q = 0.29 and q = 0.39, respec- 
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wavelength (A) 



Figure 10. Pa7 and He I 1.083 fim line profiles (down) and telluric tem- 
plate (up). Orbital (left) and supercycle (right) phases are indicated. Note 
the discrete absorption features at the blue wing of the Pa7 line. 



lively. As the ratio between the stellar radii is ~ 1.85, we get Ri 
0.16a. 



4 DISCUSSION: PROPOSED SCENARIO AND 
RELEVANCE FOR DPVS 

The previous analysis could indicate that 

OGLE05 155332-6925581 is a semi-detached binary of the 
Algol type. In these binaries the less massive star appears more 
evolved than the more massive one, a fact that is explained if the 
former was initially the more massive star, evolved first and then 
has become the less massive component by processes of mass 
transfer and mass loss. The system OGLE05155332-6925581 ap- 
parently fulfills this criteria, since it shows g < 1 and the secondary 
star in an advanced stage of evolution filling their Roche Lobe. 
This interpretation, if confirmed, could be relevant, since very few 
Algols have been found in the OGLE database of the Magellanic 
Clouds (Wilson 2004). 
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Figure 11. Pa/3 line profiles (down) and telluric template (up). Orbital (left) 
and supercycle (right) phases are indicated. Note the discrete absorption 
features at the blue wing of the Pa/3 line. 



Figure 13. Radial velocities for the discrete absorption components. We 
label weak (W) and strong (S) lines. Dashed lines indicate very weak or 
absent fines. 
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Figure 12. The Pa/3 line profile at phase 0.47 (superphase 0.22) along with 
the synthetic binary profile built considering vsini\ ■ 
= 100 km s~^ and the parameters given in the text. 



Figure 14. Normalized average SMARTS spectrum with labels and equiva- 
lent widths (in A) for the main lines. The cross indicates a bad pixel region. 



4.1 On the circumprlmary disc luminosity and mass ratio 

We have obtained additional constraints on the mass ratio us- 
ing the knowledge of the distance, orbital period, inclination 
angle, half-amplitude of RV, temperature of the secondary, as- 
sumption of spherical geometry, semidetached configuration of 
the system and calculated fractional light of the components. 
From the later and observed AIv = —3.52 we estimate Mv,i 
= —2.52 and Mv,2 = —2.96 as the upper limits for the stellar 
magnitudes. They should correspond to the stellar magnitudes in 
the case of pure stellar models and zero disc light contribution. 
Then we calculated the theoretical Mv,i{Mi{q,K2,i, P), Z) and 
Mv,2{fi2[P),T2, M2{q, K2, i)) for a set of K values assuming a 
main-sequence mass-luminosity relationship plus stellar evolution 
models (Bertelli et al. 1994) for the first case and the mean den- 
sity of a Roche-lobe filling star plus bolometric corrections given 



by Bessel et al. (1998) in the second case. We used Z = 0.008 and 
E{B-V)=0.n. 

Fig. 19 shows our results in the more stringent case of the low 
temperature limit of T2 = 12 500K. To interpret the diagram we 
assume that the third-light contribution is neglected but not those 
of the circumprlmary disc. In this case, it is evident that the g > 1 
case is rejected, as can be easily seen in the behavior of the solid 
K2 lines that cut the AIv,2 horizontal line at q values lower than 
~ 0.75, that corresponds to the K2 value of 120 km s~^. At this 
point we note that the diagram can be used as a diagnostic for the 
circumprlmary disc luminosity. For the case of K2 = 120 km s~^ 
and looking the corresponding dashed line we observe that the disc 
contribution should be about 2.5 magnitudes to the disc-l-primary 
subsystem. On the other hand, if K2 is 170 km s~^ then q < 0.4 
and the disc contributions is about 1 magnitude to the disc-l-primary 
subsystem. In all cases we observe the disc more luminous than the 
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Figure 15. Normalized high resolution MIKE spectra around Ho, H/3 and 
H7 (from up to down) taken at orbital phase 0.756. The velocity zeros rep- 
resent the laboratory rest wavelengths. Note the narrow emission feature 
inside the weakened Hq absorption and the asymmetry of H/3 and most no- 
tably H7. The vertical dotted line indicates the velocity of the Ha emission 
maximum. 



Figure 17. Long-term OGLE III /-band light curve along with differential 
spectrophotometric magnitudes mj (diamonds) and mjplus (filled circles). 
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Figure 16. Radial velocities and /-band magnitudes folded with the orbital 
period. The best sinus fits are also shown. The vertical bars indicate the sin- 
gle measurement errors derived from the RMS of the standard star radial 
velocities. Note the 7 shift between optical and infrared lines. MIKE veloc- 
ities correspond to the H7 minimum and H7 gaussian fit minimum from the 
lower to the higher velocity. The velocity of the peak of the Ha emission is 
indicated by a vertical dash. 
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Figure 18. Simultaneous BVI magnitudes obtained at maximum and si- 
multaneous JHK magnitudes relatives to the B magnitude. The magni- 
tudes have been dereddened using Ay = 0. 17. We also show solar metallic- 
ity models from Bessell et al. (1998, for log g = 3.0) and a low metallicity 
model from Bertelli et al. (1994). The bars indicate the 90% probability loci 
for JHK magnitudes at maximum based on Monte-Carlo simulations, i.e. 
the knowledge that they were obtained from the average of 10 randomly dis- 
tributed images (Ita el at. 2004) and the assumption that the OGLE II light 
curve can be used as a seed for randomly generated infrared light curves. 
The figure shows that the best model for the BVI magnitudes at maximum 
is that with T = 20 000 K, and that there is a color excess already at the 
I-band that is prominent in redder passbands. 



primary star which could explain the general absence of observed 
primary star features in the spectra. This finding cautions about the 
assumptions of our LC model method and its results. The primary 
eclipse could be mainly from the disc, not the primary star. In this 
sense, the system could be similar to /3 Lyr (Wilson 1974). We note 
that the g < 1 value derived here is independent of our LC model- 
ing and robust against the presence of the circumprimary and outer 
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Figure 19. The g-magnitude diagram showing dynamical and model 
constraints for both stellar components. The horizontal lines labeled 
Mg, Mv,2 and My^i indicate the absolute magnitude of the system, sec- 
ondary and primary star, respectively, for the case of a pure stellar binary 
(i.e. non disc). In the real world these are upper limits since they include un- 
known contributions from the circumprimary and circumbinary discs (this 
uncertainty is showed by arrows). The oblique dashed lines represent al- 
lowed primary star My for three different K2 values from Table 4. The 
oblique solid lines show the secondary star My values for the same set of 
RV half amplitudes. These data comes from modeling of the secondary star 
as a filling Roche-lobe component with fixed T2 = 12 500 K (lower limit). 
See Section 4.1 for details. 



discs. On the other hand, the possible non-neglectable contribution 
of the circumbinary disc makes the above estimated magnitudes 
only upper limits for the circumprimary disc contribution. In the 
high K2 case, the circumprimary disc could be comparable in lu- 
minosity with the primary star and with the circumbinary disc. 



4.2 Mass accretion and P 

If the disc is the main source eclipsed during main eclipse, and its 
brightness is fully accretion powered and comparable in luminosity 
to the secondary star, then the mass accretion rate should be related 
to its luminosity through: 



^M, 



(5) 



Using an approximated mass-radius relation for main sequence 
stars Ri = i?Q(Mi/MQ)"-^, the Mv = -2.83 for the secondary 
star derived above and BC = — 1.2 for a 14000 K star, we obtain 
Ml ^ 1.03 X 10"^(Mi/Mq)"°-^ Mq yr"\ weakly dependent 
on the primary star mass. For Mi > 1.8 Mq we get Mi > 8.6 
xlO"'' Mq yr^^. In the conservative case of mass transfer, we 
should observe an orbital period increase given by: 



Pq _ 3Mi(Mi-M2) 
IP^ ^ M1M2 



3Mi(l-<;) 
M2 



(6) 



(e.g. Hilditch 2001). For q ^ 0.35 and M2 < 2.8 Mq we should 

observe p > 6 xlO^^ yr^^, being the observed upper limit 7.5 
X 10^® yr^^- We propose in Section 5 an explanation for the rela- 
tive constancy of Po- 



4.3 Possible evolutionary stage 

The initially more massive star of a close binary fills their Roche 
lobe during one of three possible phases (Kippenhahn and Weigert 
1967): during core hydrogen burning (Case A), H-shell burning 
(Case B) or He-shell burning (Case C). The fact that we observe 
a hydrogen secondary star in OGLE05155332-6925581 would in- 
dicate that the system is the result of a Case A or Case B mass 
exchange. Calculations of the evolution of a close binary compo- 
nent of 9 Mq during Case A mass exchange show a short phase of 
fast mass loss (at rates ~ 10^"^ MQ/year) from the initially more 
massive star lasting ~ 6 x lO** years, after which the star loses a 
large fraction of their initial mass, reaching 3.73 Mq (Kippenhahn 
and Weigert 1967). During this phase the system loses matter and 
angular momentum through the L2 point forming a Keplerian ring 
located at ~ 2.25 x a (Soberman et al. 1997). As this rings car- 
ries a lot of mass and angular momentum the result is an important 
decrease of the orbital period. As we do not observed this rapid 
period change in OGLE05155332-6925581, the system is proba- 
bly after this rapid stage, in the g < 1 regime, consistent with our 
photometric mass ratio and conclusion of Section 4.1. 



4.4 On the long-term variability and circumbinary disc 

The non-eclipsing nature of the long variability places their ori- 
gin outside the binary system, probably in a circumbinary disc. It 
is highly likely that this disc is also responsible for the observed 
infrared excess. We note that this disc is probably not the fossil 
of the disc formed during the phase of rapid mass loss by the bi- 
nary, since we found evidence that new material is being constantly 
supplied by the binary. We propose that gas is being expelled out- 
side the binary through the L2 and Ls points. The DACs pattern of 
Fig. 13 strongly supports this view and indicates that the gas accel- 
erates through the streams in open arcs around the binary in oppo- 
site directions to the orbital motion feeding the circumbinary disc. 
It is possible that the remanents of earlier evolution stages form 
an extensive structure around the system modifying the observed 
properties of the object (e.g. colors). Further observations in the far 
infrared and radio wavelengths should be useful to check this view. 

The loop encountered in the CM diagram could provide a key 
to understand the cause for the long term variability. Similar loops 
have been found in SMC Be stars (de Wit et al. 2006), being inter- 
preted in terms of Bremsstrahlung emission arising from cycles of 
disc formation and dissipation around Be stars. De Wit et al. (2006) 
showed that loops in the CM diagram appears when optically thick 
discs are formed by stellar mass ejections and then dissipate as op- 
tically thin rings into the interstellar medium. In the same way, the 
CM loop observed in OGLE05155332-6925581 could indicate the 
formation of an outflowing optically thick circumbinary disc and 
their posterior dissipation in the form of an optically thin ring into 
the interstellar medium. The correlation found in DPVs between 
the orbital and long period (Mennickent et al. 2003) suggests that 
the mass loss is modulated by some kind of instability inside the 
binary system. 

It has been shown by Packet (1981) that during mass exchange 
in a closed binary system, enough angular momentum is transferred 
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toward the mass accreting star to spin it up to its critical rotational 
velocity after gaining only a small percentage of its original mass. 
Hence, the primary of OGLE05 155332-6925581 could be rotating 
at their critical velocity (see also de Mink et al. 2007). Our obser- 
vations of a rapidly rotating primary support this view. In this case 
additional gas supplied by the secondary could form a disc around 
the primary. The residuals observed in the LC models, the Ha emis- 
sion and the changing He I lines, could be evidence of a disc plus 
interacting region. We propose that this disc fills the Roche lobe 
of the primary and escapes through the L3 point. If the system has 
a low mass ratio, the disc could grown beyond the 3:1 resonance 
radius (Lubow 1992, Murray 1996) and precession could help to 
funnel material through the L3 point (Mennickent et al. 2003). The 
behavior of the helium lines, which are detected primarily during 
long cycle minimum shows that the interacting hot region or its 
visibility is affected during episodes of strong mass loss. 



5 CONCLUSIONS 

In this paper we have analyzed the eclipsing double periodic 
variable OGLE05155332-6925581 arguing that it could be a 
clue to understand the phenomenon of Double Periodic Vari- 
ables. We find a set of physical parameters best fitting the mul- 
tiwavelength light curves and find a much improved value of 
the orbital period. Our LC and dynamical analyses indicate that 
OGLE05 155332-6925581 is a semi-detached binary with inter- 
mediate mass components where the less massive star transfers 
matter onto the more massive star. We detected evidence for a 
luminous disc of gas around the primary as in /3 Lyr. However, 
OGLE05 155332-6925581 does not show changes in the orbital 
period as P Lyr does. In our view the system is found after a rapid 
mass exchange associated to Case A or Case B mass transfer, where 
the secondary star appears undermassive and overluminous due to 
the erosion of their outer layers by the mass loss process. 

We find evidence for a variable circumbinary disc, that should 
be the cause of the long-term periodicity and the infrared ex- 
cess. Remarkably, the long-term variability follows a loop in the 
CM diagram (Fig. 5) indicating quasi-cyclic episodes of disc cre- 
ation/dissipation; the mass loss being modulated by some still un- 
known process in the binary. The remarkable discovery of DAC 
RVs following a saw-teeth pattern during the orbital cycle strongly 
supports the view that mass is being lost from the system through 
accelerated streams arising from the Lagrangian L2 and I/3 points 
(Fig. 13). This should be the source of mass for the circumbinary 
disc. The constancy of the orbital period could indicate that the op- 
posite effects of period increase and period decrease due to mass 
exchange and mass loss cancel in a long time scale. The supercycle 
clock could act like a safety valve, allowing the mass to cumulate 
around the primary during a supercycle and then be expelled away 
from the binary at the end of the supercycle, with almost no long- 
term effect on the system orbital period. 

We conclude that OGLE05 155332-6925581, and probably all 
DPVs, are intermediate mass semidetached Algol type binaries in a 
evolutionary stage characterized by mass exchange, mass loss and 
circumbinary discs. DPVs looks more massive, bluer and more lu- 
minous than classical Algols; their study is potentially important to 
constraint non-conservative models of binary star evolution. 

Despite of the significant advances in the knowledge of 
DPVs provided by the study presented in this paper, we are 
certain that problems with the complete interpretation of the 
phenomenon still persist. The complex and variable nature of 



OGLE05 155332-692558 1 requires better quality and more exten- 
sive observational material than we have collected. Firstly, we do 
not know what is the mechanism driving and regulating the long 
cycle variability. In addition, the identification of detected spectral 
features with specific places in the system and its environment is 
still uncertain, especially for the Ha emission. We do not know 
how important is the role of the mass transfer and fast rotation of 
the primary star in the mechanism driving the long cycle variability. 
The interesting 7 shift observed between optical and infrared lines 
of the secondary still has no explanation. The stellar parameters for 
the primary obtained from our LC modeling are probably biased 
by the presence of a poorly known circumprimary disc. Finally, 
the basic properties (mass, size, temperature) of the circumbinary 
disc are still unknown. The recent discovery of bright DPVs in the 
Galaxy (Kolaczkowski & Mennickent, in prep.) open new oppor- 
tunities to our research. We are going to extend our project for in- 
tensive, multi-instrumental study of these very promising targets. 
It will permit to reach the level of accuracy and time resolution 
essential for constructing a comprehensive model of their activity. 
On the other hand, the Magellanic Cloud DPV sample is 
unique in their completeness and homogeneity. For these systems 
we roughly know distance, total absolute magnitude and luminos- 
ity; we will continue working on this sample to provide a sound 
basis of comparison for the Galactic sample with the goal of get- 
ting insights on the influence of the metallicity on the DPV prop- 
erties. We plan also make a detailed comparison with the Galactic 
classical Algols to enlighten the DPV-Algol connection. 
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